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@ Das Verfahren nutzt die Gradienteninformation zur Mi- 
nimierung gestorter Zielfunktionen. Dabei wird der Gra- 
dient oh ne zusatzliche Zielfunktionsauswertungen mittels 
RLS-AJgorithmus geschatzt und mit einem stochasti- 
schen Mmtmierungsverfahren kombiniert. 
Anwendung des Verfahrens ist die Optimierung von Mo- 
del I para metern eines Kuhlstrecken model Is zur Minirnie- 
rung des Modellfehlers. 
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Beschreibung 




In der Industrie werden haufig ProzeBmodelle zur Beschreibung tcchnischcr Anlagen cingeseizt. Dabei and noch un- 
bekannte ProzeBparameter zu bestirnmen und die Werte ungenau bekannler Parameter sind zu verbessern. 
5 Zur Ldsung iiberbestimniter Gleicbungssysteme, die keine exakte L6sung besitzen, kann vorteilhaft mitlels globaler 
stochastischer Optiimenmgsverfahren eine geeignete Ldsung numerisch durch einen Rechner ermittelt werden. Dies ist 
vor allem bei ProzeBmodellen sinnvoll, die eine groBe Anzahl Parameter besitzen. Standardalgorithmen zur Minimie- 
rung des Modellfehlers scheitern an zahlreicben lokalen Minima, an Dumerischen Problemen bei der Gradientenberech- 
nung oder an der Rechenzeit 
10 Die Aufgabe des erfindungsgemaBen Verfahrens besteht darin, ProzeBparameter numerisch zu optimieren. 
Diese Aufgabe wird durch das \ferfahren des Patentanspruchs 1 gelost 

Zur ErrnMung oder Verbesserung dieser ProzeBparameter werden an der techniscbe Anlage zahlreiche MeBdaten auf- 
genornmen und die AusgangsgroBen der Anlage mit den Modellausgaben vergiichen. Durch Minimierung des Modell- 
fehlers erhalt man die gesuchten Parameterwerte. 
15 Das erfodungsgemaBe Verfahren kombiniert ein stochastisches Optimierungsverfahren mit einem in der Regelungs- 
technik und Nachrichtenabertragung bekannten Schatzalgorithmus, der zur Schatzung einer Abstiegsrichtung (negativer 
Gradient) umformuliert wurde. Das stochastiscbe Optimierungsverfahren setzt dabei voraus, daB die Zielfunktion zwei- 
mal stetig differenzierbar ist 

Da in der Praxis diese Voraussetzung an die Zielfunktion selten erfUllt ist, ist die numerische Gradientenberechnung 
20 oft problematisch. Dies triflt vor allem auf Optimieiungsprobleme zu, bei denen die Zielfunktion beispielsweise durch 
das Softwarepaket des Kunden spezifiziert ist Bei der dem erfindungsgemaBen Verfahren zugrunde liegenden Minimie- 
rung geniigt die Kenntnis eines Gradienten eines zweimal stetig differenzierbaren Anted Is. Dieser Gradient kann mit dem 
RLS-Verfahren robuster als mit vergleichbaren Methoden der numerischen Mathematik geschatzt werden. Dabei ist es 
vorteilhaft, daB der Gradient ausschlieBlich aus zuruckliegenden Funktionsaufrufen ermittelt wird. 
25 Es ist vorteilhaft, wenn der Schatzalgorithmus zur Schatzung der Abstiegsrichtung mit einer Zuf allssuche kombiniert 
wird, da so die Suche nach dem globalen Minimum robust wird gegen Lokale Minima. 

Auch ist es vorteilhaft, wenn ein Vergessensf aktor berucksichtigt wird, um den EinfluB weit zunlckliegender schlech- 
ter Werte bei der Iteration uber alle ermittelten Werte abzuschwachen. 

Weiterbildungen des erfindungsgemaBen Verfahrens ergeben sich aus den abhangigen Anspriichen. 
30 Das erfindungsgemaBe Verfahren wild anhand der folgenden Figuren weiter erlautert 
Es zeigen 

Fig, 1 ein Blockdiagramm, das die einzelnen Komponenten des Verfahrens darstellt 

Fig. 2 eine Skizze, die die Funktion des erfindungsgemaBen \ferfahrens veranschaulicht, 

Fig. 3 eine Skizze, die die Funktion der Zufallskomponente veranschaulicht, 
35 Fig. 4 ein Blockdiagramm, das eine Anwendung des Verfahrens darstellt. 

Fig. 1 zeigt das der Erfindung zugrundeliegende Prinzip. Das dargestellte Optimierungsverfahren (Block la) erzeugt 
einen Wert x e R° (x-Wert XW), an dem eine Zjelfunktionsauswertung vorgenommen wird (siehe Block lb). Es ist der 
Gradient an dieser SteLLe x erforderlich, um einen neuen x-VfertXW zur Funktionsauswertung zu bestimmen. Diese Be- 
stimmung eines x-Wertes XW kann iterativ erfolgen, bis der k-Wert XW mit dem kleinsten Funktionswert gefunden ist 
40 Der Gradient stent i.a. nicht zur Verfiigung, so daB er numerisch berechnet werden soil. Bei dem erfindungsgemaBen \fer- 
fahren wird der Gradient mittels zuruckliegender Funktionsauswertungen rekursiv geschatzt (sie Block lc), Dieser Gra- 
dient bestimmt die Suchrichtung SR nach dem globalen Minimum in der nachsten Iteration des globalen stochastischen 
Optimierungsverf abrens (la). 

Dabei werden die x- Werte XW nach einem Zuf allsprinzip erzeugt, was eine robuste Gradientenschatzung durch den 
45 verwendeten RLS-Algorithmus erlaubt. Beispielsweise kann als Eingangssignal fur den RLS-Algorithmus weiBes Rau- 
schen verwendet werden. 

In der Praxis ist die Voraussetzung der zweimal stetigen Differenzierbarkeit der Zielfunktion selten erfuilt. Bei der glo- 
balen Minimierung genugt die Kenntnis des Gradienten des zweimal stetig differenzierbaren Anteils, der mit dem RLS- 
Algorithmus robust bestimmt werden kann. 
50 Somit sind Vbrteile des erftnduragsgernaBen Verfahrens, daB keine zusatzlichen Zielfunktionsauswertungen erforder- 
lich sind und daB die Gradientenschatzung gegen Storungen unempftndlich ist. 

An Fig. 2 soli im Folgenden das Verfahren zur globalen Optimierung naher erlautert werden. 

Als Voraussetzung sei f: R n — * R eine zweimal stetig differenzierbare Zielfunktion, die ein globales Minimum besitzt 
Die Aufgabenstellung besteht in der Bestimmung dieses globalen Minimums. 
55 Das erfindungsgemaBe Verfahren stellt eine stochastische Methode dar, die auf der numerischen Integration folgender 
stochastischer Integralgleichung basiert (siehe [2]): 

X^CD, t) = Xq - J Vf^X^ ((D, T))dT + £B(o, t) (1) , 

0 

mit a> e Q 
wobei 

65 Qeine Menge aller Ergebnisse to, 
t einen Zeitpunkt t e Rq\ 
x 0 einen Startpunkt x 0 e R°, 
X eine Abb. X: QxR^ — R°, 
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B eine Abb. B: £2xRo* R n i 
e einen vorgebbaren Parameter i 
bczeichncn. 

Dabei liegt der Wahischeinlichkeitsrauni (Q, S, P) zugrunde und die Abb. B wird als n-dimensionale Brownscbe Be- 
wegung bezeichneC Der Gradient der Zielfunktion wird durch S 



10 

und die Zufallssuche durch B dargestellL 

Die stochastiscbe Integralgleichung (1) kombimeit somit die lokale Mininrierung (Kurve des steilsten Abstiegs) mil 
der Zufallssuche (Brownsche Bewegung). 

In Fig* 2 wird die numerische Losung der Gleicbung (1) fur einen Pfad des stocbastischen Prozesses veranschaulicht 
Dabei bedeutet der stochastische Prozefl eine Sammlung 15 



20 

von D-diinensiarialen Zufallsvariablen 
X^-.Sl -» R n ,-CD W> X^(o,E) 

25 

fur jedes feste t e Ro + 

und ein Pfad des stochastischen Prozesses die Funktion 



XQ b ' fc - - ' - "XqV / 30 

fur jedes feste <£> € SI 

Die Zufallssuche ist in Fig. 2 duich einen Kreis ZSK und die Gradientenrichtung durch einen Pfeil GPf symbolisiert 
Beftndet man sich im Punkt 



so geht man zur approximativen Berechnung von 



x xq(®' t + T ) mit T ^ 0 sr eTnas " TVf ( x xh(^' ^)) 



35 



40 



in Richtung des negaliven Gradienten, verrauscht dann die Gradientenrichtung gemaB 45 

s(b(cD, t + T) - B((D , t)) , 
erhak den neuen Punkt 50 

und setzt die Iteration fort 55 

Die Wahl des Parameters £ wird vorgebbarnach folgenden Kriterien vorgenommen (siehe [2]): 
Zeigt sich bei der Betrachtung der bisher berechneten Iterationspunkte, daB bei der numerischen Losung von (1) die 
Kurve des steilsten Abstiegs 



60 



fur ein festes ©eft dominiert, so ist £ zu erhohen . Ist hingegen bei der Betrachtung der bisherigen Iteratioaspunkte die 65 
Zufallssuche 
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5 fiirein festes <D e Q dominant, so ist e zu veningem. 

Somit ist der Parameter e zur Steuerung fur cine Regularisierung zwischcn lokaler Minimierung und Zufallssuche zu 
vcrwcnden (Regularisierungsfaktor). 

Zur Losung von Gleichung (1) konnen numerische \ferf ahren fur Anfangswertprobleme gewohnlicher Differential- 
gbachungen eingesetzt werden. Da aber der ALgorithmus fur Zielfunktionen ausgelegt ist, bei denen weder der exakte 
io Gradient noch die exakte Hessematrix zur Verfugung stehen, wurde als Verf ahren zur lokalen Minimierung das explizite 
Eulerverfahren (Gradientenverfahren) gewahlt Dabei ist wesentlich, daB der benotigte Gradient mittels ULS-Algorith- 
mus geschatzt wird. 

Fur die Gradientenschalzung wird von folgender Gleichung ausgegangen: 

mit y = f (x) e R 



wobei 

20 y einen Funktionswert von x e R n , y 0 einen Schnittpunkt der Ebene mit der y-Achse und 
g die zu schatzende Suchrichtung 
bezeichnen. 

Im Folgenden werden k+1 Punkte y<> y lP y 2 , . . yk betrachteL Diese k+1 Punkte konnen beispielsweise durch Itera- 
tion bei einer numerischen Optimierungsaufgabe erzeugt worden sein. Dabei bezeichnen k die diskrete Zeit, den ak- 
75 tuellenPunktundx^!,. . ^xj^XoweiterzuriickliegendeStutzpuntornitXi e R°,i = 1,2,. . k. Durch Einsetzen dieser 
Punkte in (2) ergibt sich das Gleichungssystem 
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yi 



i 



(3) , 



35 wobei gt e R° die anhand der genannten k+1 Stiitzstellen zu schatzende Suchrichtung und y^ e R n den dazugehori- 
gen Schnittpunkt der Ebene mit der y- Achse bezeichnen. Eine Subtraktion der einzelnen Gleichungen in (3) bewirkt cine 
Elimination der Konstanten y^, die fiir die Gradientenschatzung nicht von Interesse ist. 



( 



yi - yo 
yt. - yi 



x'2-x'! 
*k 



3k 



(4) 



VYk - yic-i/ 

' dk ' 

Der Gradient g k stellt dabei den zu schatzenden Gradienten V f(Xfc) der Zielfunktion f an der Stelle x k dar, Zum Zeit- 
punkt k+1 lautet die Gleichung (4) dann: 



f yi - yo ^ 
y2 - yi 

Vyk+l ~ Yk/ 
d k+l 



x'2-x'i 
^"k+i-x-jj 

v 1 ■ ■ 

*k+l 



9k+l mit 3k+l e * 



n 



(4a) 



Das Gleichungssystem (4a) enthalt i.a. mehr Stutzpunkte (Gleichungen) als Parameter, also gilt k+1 > n. Wie ein sol- 
ches uberbestimmtes Gleichungssystem in dem erfindungsgemafien \ferfahien gelost wild, wird im folgenden Abschnitl 
erlautert 

Der LS-Algorithmus (Least-Square- Algorithmus) wird dazu verwendet, uberbestimmte lineare Gleichungssysteme 
veranderlicher Gr6Be zu losen. 

Fiir ein festes k aus der Menge der naturlichen Zahlen N betrachtet man das Gleichungssystem 
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M kPk = *>k 
mit 



(5) 



Mk = 



m« 2 



e H k ' n , 



mi 6 R n / 



i=l, 2, . . ,k, 



e R und 



bk = 



02 



Pi e R; k > n . 



In der voriiegenden Darstellung werden alle \fektoren als Spalteovektoien dargestellt 
I.d.R existiert kein p^, das die Gleichung (5) exakl erfiillL Das Residuum 

r (Plc) : = V" M kPk 

ist also fur kein zum Verschwinden zu bringeix Als Ersatz fur die exakte Losung wird ein gesucht, das das Resi- 
duum moglichst klein werden LaBt, z. B. im Sinne der euklidischen Norm soil ein p k so bestimmt werden, daB gilt 



Pk)I 2 * WPkJL V Pk e R 1 
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20 
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Dies wird als das lineare Ausgleichsproblem oder die Ausgleichung der Widerspriiche nach der Methode der klcinsten 
Quadrate bezeichnet 

Weiterhin kann derEinfluB writer zuruckltegender Daten verringert werden, damit der zu schatzende Parametervektor 
Pk den vorgegebenen Parametero des Gleichungssystems folgen kann. Dieses "Vergessen" bewirkt, daB die Daten mit ei- 35 
nem um so klcineren Gewicht versehen werden, je weitcr sie zuruckliegen. 

Ein Mogb'chkeit hierfur ist das multiplikative \fergessen. Mittels der Gewichtungsmatrix 



ft 



Dk - 



d2 



mit d k = 1, d k . v = <* v , v =1,2, . . ,k-l 



40 



45 



wird die mittlere quadratische Fehlerfunktion gebildet: 

EkCPk) - = ^(^-^Pk)'Pk(bk-MkRc) (6), 50 

wobei a ein fester vorgebbarer Parameter im Interval! 0 < a ^ 1 ist und dazu verwendet wird, die einzelnen Summanden 
in Efc(p k ) unterschiedlich zu bewerten, und zwar um so schwacher, je weiter ihr Argument k-v gegenuber k zuruckliegt 
Mit dem Parameter a kann also erreicht werden, daB die zuruckliegen den Summanden mit zunehmender Zeit k immer 
mehr "in Vergessenheit" geraten. Dieser Effekt tritt nicht auf, wenn a = 1, & h. Dk =I k V k e N gewahlt wird; dies stellt 55 
den Spezialfall der ungewichteten quadratischen Fehlerfunktion dan 

Es soil nun pt so bestimmt werden, dafl E k (p k ) minimal wird. Dazu wird der Gradient der Fehlerfunktion gleich Null 
gesetzt: 

VE k(Pk) = 0 ~ M' k DfcMkPk-M'k D k b k = 0. 
Somit gilt, falls \T k DkM k invertierbar ist: 

Pk = R*f 1 'k (V 
mit 
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Rk = M' k D k M k (8), 
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10 



15 



20 



r k = M , k D k b k (9). 

Setzt man in den Gleichungen (5), (7), (8) und (9) = pj, = g k und ^ = 4, so ist das Problem der Giadienten- 
schitzung aus GLeichung (4) prinzipiell gelost Allerdings ist die durcfa die Gleichungen (7), (8) und (9) beschriebene Be- 
rechnung noch nicbt brauchbar, da bei Hinzunahme eines weiteren Vektors m' k+1 und eines weiteren Skalars 0^ erneut 
eine Matrixinversion mit Berechnung von R^i anfallt Die GroBe des Gleichungssystems wachst mit steigendem k aa 
Der nachfolgende Algorithmic umgeht diese Schwierigkeit nrittels rekursiver Berechnung. 

Bei dem RLS-Algorithmus wird das urn eine weitere Gleichung erganzte lineare Gleichungs system 

M k+1 Pk+1 = kk+l 
rait 

M k +1 = {J^ J i Mk ^ R k ' n , t^+l e R n und 
bk+i - (p^J i bk e R n , Pk +1 e R 



betrachtet Dabei stellt pk +1 den neuen zur Zeit k+1 gesuchten Vektor dar, mit dem jetzt k+1 Bedingungen bestmpglich zu 
erfulleo sind DLe Zeilenzahl der Matrix M k + 1 und die Lange des Vektors bk +1 sind also gegenuber M k und bt urn " 1" an- 
gewachsen. 

25 Das wesentliche Merkmal des RLS-Algorithmus ist nun die rekursive Ermittlung der inversen Matrix RjJ, . Daher ist 
das Ziel der nacbfolgenden Betrachtungen, eine Losung folgender Form zu finden: 



( neuer 
30 Vschat zwertJ 



^ ( alter ^ ( Korrektur f Korrektur -A 
:) " Vj5ch&tzwertJ Vvektor J Uaktor J 



In den foigenden Schritten wird das Aufstellen der rekursiven Gleichungen naher erlautert. Aus den Gleichungen (7), 
(8) und (9) folgt fur den Zeitpunkt k+1 

35 Pk+1 = *k+l * *+l {10) ' 

R k+1 = Mk+l^k+lMk+l (11), 
2fc +1 = Mk+iDfc+xbfc+x (12). 

40 

Durch Ausmultiplizieren ergeben sich aus (11) und (12) folgende Gleichungen: 
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k+l k+1 _ y 
R k+1 = X a ~ V mvm' v 
v=l 

mit m' v e R (v-te Zeile von Mfc+i) 
k 

o R k+1 = a 2 a k ^vm'v+ttk+im'k+i 
v=l 

o ^k+i = aRfc + mjc+i^'k+x <13i 



und 



Gleichung (13) stellt eine Recursion fur die Berechnung der Matrix R^ +1 dar. Daraus MBt sich unmittelbar eine Recur- 
sion fur die inverse Matrix R*» 1 , bestinunen (siebe [1]): 



R "l _ 1 

*k+l 



-1 1 -i> 

i _ Rk • "k+i • m k+ i • 



(15) , 



a + "hc+l " R k * "k+l 

erhait man 

Kfc+l - ^( R k 1 - h k+l ^c+i Rfc" 1 ) <16) . 



to 



15 



20 



k+l ^ 
v=l 

mit JJ V € Rn (v-te Konrponente von bfc +1 ) 

k 25 
<=> *k+l = ct 2 a ^v^v + Pk+l m k+l 
v=i 

<=> *k+l = + ftc+l m k+l (14) . 



30 



35 



d. h. man kann ohne Matrixin version aus der inversen Matrix R*~ l dcs vorherigen Schrittes berechnen. 40 
Mit der Abkurzungen 

h k+l - , 1= 



45 



50 



Fiir den zu schatzenden Parametervektor ergibt sich durch Einsetzen von (14) und (16) in (10) und nach weiterenUm- 
formungen die Gleichung (siehe auch [1]) 55 

Pk+l = Pk 4 n k+i(0k+t- m k+l Pk) 
bzw.mite k = p k+1 -m' k+l p k 

Pk+i=Pk-*"hk+i^ (17). 



60 



Fiir das rekursi ve Verfahren werden vorteilhaft die Startwerte pj und Rj -1 benStigt. Als Startwerte werden diejenigen 
Werte verwendet, die einen moglichst geringen EinfluB auf das Schatzergebnis nehmen: Es empfiehlt sich als Startwerte 
Pj » OundR] -1 mitTi > 1 zu verwenden. Die dadurch entstebende Ungenauigkeit beeinffuBt den Vektorp^. Inder 65 
Praxis ist es vorteilhaft, den Parameter T\ moglichst groB zu wahlen. 

Der RLS-Algoritbinus gewinnt seine Information iiber die Datenmatrix M^ und den Datenvektar und speichert 
diese Informadon in der Matrix R^ 1 . FUhrt man eine Eigenwert-Eigen vektor-Zerlegung der Matrix R k durch, so ergibt 
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sich 

Rk = VMc^k 
s mit Aj^ = diag^X^, ^2^' * • * ^nk) ' 

wobei Xi k ; i = 1 T 2, . . n die Eigenwerte der Matrix H± darstellen und die Matrix 

aus den Eigenvektoien vt: i = 1, 2 n der Matrix Rfc besteht 

Aus den Gleichnngeo (13) und (18) folgt fur deo Fall des multiplikativen \fergessens 



(18) 



(19) 



15 



20 



25 



30 



(aX 



aX 



■P k +m }c+1 m , k + i 



mit 



(aX 



Ik 



Tk 



(20) 



Fiir multiplikalives Veigessen gilt also folgende Bezichung fur die Eigenwerte \ von Rfc und die Eigenwerte 5tj k 
R k : = ctR k 

35 = a • mit i=l,2,..,n. 

Da die Matrix R^per Definition symmetrisch und positiv (senri-)definit ist, folgt 
40 Xj k > 0 fur i = 1, 2, . . ., n 
und 

T'kT^I,,. 

Durch Inversion von erhalt man 

( i 
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aX 



'Ik 



aX 



'2k 



Tk 



(21) 



60 bzw. fiir die Eigenwerte 
1 = 1 
^±k aXi k 



fur i=l, 2, . . ,n 



(22) 



65 



Es ist vorteilhaft, wenn anstelle der Gleichung (22) folgende Gleichung verwendet wird: 
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^ X *k< 



fur i=l, 2, . . ,n 



(23), 



wobei je nacb Wahl von 8 verschiedcne Einstellungen filr das Vergessen moglich rind: 
- mulliplikatives "Verg 3 sen wie in Gleichung (22) fiir 5 = 0, 

• additives Verges sen tni|. 5=1: ^A- 1 



v *k 



l - a 
a 



10 



Durch Einsetzen der Diagonalelemente in Gleichung (20) mit den Kehrwerten der Eigenwerte gemSB Gleichung (23) is 
ergdbt rich 



a 



a <4c 



Tk ( 24 > • 



1 + — "4 7 
a 

Durch Ersetzen von anstelle von CtRj- in (13) ergibt sich 

R k+1 = fi k + m k+l m, k+l • 

Ersetzt man die Diagonalelemente der Matrix aus (21) durch die Eigenwerte der Gleichung (23), so eigibt sich 



f 



a 



V 



ik 

l | l-a l 



(25) . 



-1 



Aus [1] folgteine neue lekursive Beziehung 

p-1 _ p-1 _ gfc 1 ' ™k-H * m 'k+i ' 
K k+1 ~ K k , ■ ^-l 

1 + tttfc.^ • Kfc • nik +1 

Mit der Abkiirzung 



(26) . 



*-l 



h k+l = 



k ' Mk+l 



1 + m, 



^k 1 ' ^+1 



*k+l 

wird die Rekursion berechnet: 



(27) 



(28) . 
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chten GroSe p^ +1 setzt man wie beim multiplikat^TO \ 



Zur Berecfanung der^Hchten GroSe p^ +1 setzt man wie beim multiplikaflTO Vfergessen Gieichung (17) mit der Ab- 
kiirzung 

e k =p k+1 -m' k+l p k (29) 
zu 



Pk+i = Pk + h k+i (30). 

10 Nachfoigend wild nochmals zusammengefafit, wie mittels RLS-Algorithmus der Gradient einer Zielfunktion f an ei- 
ner bestimmten Stelle Xj bestimml wird. , t 

Die Iteration starlet im Schritt k= l.Es werden die Eigenwerte und Eigenvektoren der inversen Matrix bestimmt 
und die Transformationsmatrix T k aufgestellt Dann wird nach Gieichung <25) die inverse Matrix R k berechnet Die 
Parameter a und 8 sind vorgebbar gemaB Gieichung (23). Daraufhin sind die Abkurzungen und h k+1 zu berechnen 

15 nach den Gleichungen (29) und (27). Dabei sind nV k+1 durch x' k+ r* k Pk+l durch Vk+l^k g^maB Gieichung (4a) zu 
ersetzen. Danach wird p^ +1 laut Gieichung (30) mittels pt und den Abkurzungen ek und h k+t ausgewertet SchlieBhch 
wird die inverse Matrix R,^ gemaB Gieichung (28) bestimmt und gegebenenfalls mit dem nachsten Iterationsschritt fort- 
gefahren. 

Als Ergebnis stellt p^ den gesuchten Gradienten g k+x zur Zeit k+1 gemaB Gieichung (4a) dar. 

20 In Fig. 3 ist eine Skizze dargestellt, die das erfindungsgemaBe Verfahren anhand einer Zielfunktion f(x) mit lokalen 
Minima zeigt. Durch die Kombination mit der Zufallssuche endet das Verfahren nicht in einern lokalen Minimum LM, 
sondem springt unter Beriicksichtigung der Brownschen Bewegung und des negativen Gradienten aus dem lokalen Mi- 
nimum LM in Richtung globales Minimum GM. Das erfindungsgemSSe Verfahren durchiauft also im steilen Bereich die 
lokalen Minima und bleibt im flachen Bereich Qber mehrere Iterationen hinweg in dem globalen Minimum. Der Pf ad PF 

25 deutet die Suchrichtung des erfindungsgemaBen ^rfahrens an. 

Fig, 4 zeigt ein Blockdiagramm, das ein Anwendungsbeispiel fiir eine globale Optimierung mit Gradientenschatzung 
darstellt. Als Anwendungsbeispiel dient ein Kuhlstreckenmodell KSM der Stahlindustrie. Auf der Basis von MeBdaten 
MD werden Modellparameter optimiert, d. h. der Modellfehler MF unter Beriicksichtigung der Gesamtheit aller MeBda- 
ten MD minimiert 

30 Dazu werden die MeBdaten MD in das Kiihlstreckenmodell KSM eingegeben. Die vom Kuhlstreckenmodell KSM ge- 
lieferten Ergebnisse werden mit den tatsachlichen Ergebnissen der Ktihlstrecke GT verglichen und in Form eines Mo- 
dellfehlers MF der globalen Optimierung mit Gradientenschatzung GLOPGS zugefiihrt. Hier werden die Parameter des 
Kuhlstreckenmodells KSM zur Minimierung des Modellfehlers MF optimiert 
Im Rahmen dieses Dokuments wurden folgende Veroffentlichungen zitiert: 
35 [1] Unbehauen, R.: Systemtheorie. R. Oidenbourg Verlag, Wien 1993. S. 537-545. 

[2] Schaffler, S.: Unconstructed Global Optimization Using Stochastic Integral Equations, Optimization 1995, \bL35, 
S. 43-55. 



Patentanspriiche 

1 . Verfahren zur Optimierung von ProzeBparametern auf der Basis einer Gradientenschatzung auf einern Rechner, 

a) bei dem in einern Schritt k+1 ein neuer x-Wert x k+1 in Abhangigkeit von einer Suchrichtung -g k berechnet 
wird gemaB 

wobei 

fCx fc ) den zugehorigen Funktionswert aus einer vorgegebenen Funktion f bezeichnet, 

b) bei dem ein neuer Funktionswert y k+l berechnet wird gemaB 



50 



yk+1 - f ( x k+l)> 



c) bei dem ein neuer Gradient g k+1 basierend auf einern RLS-Algorithmus berechnet wird, 

d) bei dem der neue x-Wert nach einern Zufallsprinzip erzeugt wird (Rauschen). 

55 2. Verfahren nach Anspruch 1, bei dem die Schritte la) bis lc) zu mindestens einem nachsten diskreten Zeitpunkt k 

iteriert werden. 

3. Verfahren nach einem der obigen Anspruche, bei dem bei der Berechnung des Gradienten ein vorgebbarer Ver- 
gessensfaktor a verwendet wird, der x-Werte zu verschiedenen Zeitpunkten k unterschiedlich bewertet 



60 Hierzu 3 Seite(n) Zeichnungen 
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